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Oxygen Abundance Measurements of SHIELD Galcixies 

Nathalie C. Haurberg^’^’^, John J. Salzer^’^, John M. Cannon^’^, &: Melissa V. Marshall^ 

ABSTRACT 

We have derived oxygen abundances for 8 galaxies from the Survey for H I in ex¬ 
tremely low-mass dwarfs (SHIELD). The SHIELD survey is an ongoing study of very 
low-mass galaxies, with Mhi between 10®'® and 10^'® Mq, that were detected by the 
Arecibo Legacy East ALEA (ALEALEA) survey. Ha images from the WIYN 3.5m tele¬ 
scope show that these 8 SHIELD galaxies each possess one or two active star-forming 
regions which were targeted with long-slit spectral observations using the Mayall 4m 
telescope at KPNO. We obtained a direct measurement of the electron temperature by 
detection of the weak [O III] A4363 line in 2 of the H II regions. Oxygen abundances for 
the other H II regions were estimated using a strong-line method. When the SHIELD 
galaxies are plotted on a R-band luminosity-metallicity diagram they appear to sug¬ 
gest a slightly shallower slope to the relationship than normally seen. However, that 
offset is systematically reduced when the near-infrared luminosity is used instead. This 
indicates a different mass-to-light ratio for the galaxies in this sample and we suggest 
this may be indicative of differing star-formation histories in the lowest luminosity and 
surface brightness dwarf irregulars. 

Subject headings: galaxies: abundances — galaxies: dwarf — galaxies: evolution — 
galaxies: star formation 


1. INTRODUCTION 

Chemical abundance studies of the interstellar medium in galaxies allow a glimpse into the star- 
formation and chemical enrichment history of these galaxies. The chemical evolution of local- 
universe galaxies with very low-metallicity gas is particularly important as these galaxies may be 
objects of cosmological significance. These metal-poor systems represent potential local-universe 
analogs to galaxies in the early universe and can provide a better understanding of the processes of 
star-formation and gas-enrichment in the early universe. Additionally, in order for very metal-poor 
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galaxies to exist in the local universe, processes of either enriched gas outflow or pristine gas inflow 
are usually invoked, thus these galaxies can provide important constraints to understanding the 
role of gas inflows and/or outflows on galaxian evolution. 


Many authors have searched for very metal-poor or “extremely rnetal-deficient” fXMD; 12 + 


logfO/H) < 7.65) galaxies in the local universe (e.g.. Kniazev et al. 

2003: 

Ugrvumov et al. 

2003: 

Melbourne et al. 2004: Lee et a 

[. 2004 

Izotov et al 

2006: Brown et al. 

2008: 

Paoaderos et al. 

2008; 

Guseva et al. 2011: Berg et al. 

2OI2I: 

zotov et al. 

2OI2I). Some authors have extended the 

search 

these 

out to intermediate redshift (F 

iakazu et al. 2007: 

Hu et al. 2009: Xia et al. 

2012). Most of 


studies have used emission-line surveys, seeking metal-poor galaxies from large samples of strong 
emission-line sources. While these searches have led to a signihcant increase in the number of 
known metal-poor galaxies, there have been relatively few discovered in the range qualifying as 
XMD. The seeming dearth of XMD galaxies has raised the question of whether they are intrin¬ 
sically rare objects in the local universe or if they exist in large numbers but are mostly missed 
because emission-line surveys are biased toward systems with bright, high surface brightness H II 
regions (i.e., systems that have recently undergone a very strong starburst) and are relatively 
insensitive to isolated low-mass, relatively quiescent star-forming galaxies. 


A promising alternative method for finding XMD galaxies has been to study very low-luminosity 
dwarfs since there is an historically we h-established correla t ion between metallicity and luminos- 


itv among low-redshift dwarf galaxies 

(Leoueux et al. 1979: 

Skillman et al. 1989a: 

Pilvugin 2001: 

Melbourne X Salzeill2002: Tremonti et al.ll2004: Salzer et al. 1 

1005: Lee et al.ll2006: van Zee X Havnes 

20061). This method has alreadv led to the discoverv of several XMD galaxies 

Skillman et al. 

1989allbl: Ivan Zee 2000: Pustilnik et al. 

2005. 2011: Berg et a 

.. 2012: Skillman et al. 

2 OI 3 I). Despite 


the fact that there should be large numbers of low-luminosity galaxies according to the galaxy 
luminosity function, searching for these systems is very difficult as they tend to be very low surface 
brightness and difficult to detect optically. However, blind H I surveys do not suffer optical selec¬ 
tion biases and thus can uncover low-mass H I sources that potentially have low-surface brightness 
optical counterparts that may be missed in optical studies. 


It is particularly important to identify and study low surface brightness XMD sources because 
our current understanding of XMD galaxies is derived almost exclusively from blue compact dwarf 
galaxies (BCDs). However, normal low-luminosity star-forming galaxies should outnumber BCDs 
by at least an order of magnitude and thu s should be more representativ e of the XMD population. 
The Arecibo Dual Beam Survey (ADBS; iRosenberg fc Schneiden 1200(1 ) detected many low-mass 
sources with low-surface brightness dwarf irregular op tical counterpar t s. Tw elve of these ADBS 
galaxies were selected for spectroscopic follow-up by iHaurberg et al.l (j2013l l to determine their 
chemical abundances. While none of these galaxies were found to be XMD galaxies, the sample 
was generally low-metallicity and showed that HI- surveys can efficiently uncover metal-poor and 
potenti al XMD galaxies. The more recent Arec ibo Legacy Fast ALFA (ALFALFA) blind H I- 
survey (jGiovanelli et al.ll2nn5l : iHavnes et al.ll2nill ) has uncovered hundreds of low-mass H I sources, 
many with very low-surface brightness optical counterparts. The ALFALFA catalog provides a 
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promising target list of potential XMD galaxies. 


The use of H I catalogs to search for low-luminosity XMD systems has already proven successful 
with the recent discovery of the Leo P dwarf. Leo P was selected as a c andidate XMD from 
the ALFALFA catalog due to its low derived H Imass (iGiovanelli et al.ll2ni3l ) and very-low stellar 


mass inferred from optical follow-up imaging ( Rhode et al. 20131 ). Optical spectroscopy revealed 


an 


extremely low oxygen abundance of 12 -|- log(0/H) = 7.17 ±0.04 (jSkillman et al.ll2013l ) confirming 
it as an XMD galaxy and one of the most metal-poor galaxies known. 


In this paper, we report the findings from a program dedicated to determining and analyzing the 
nebular abundances of a set of low-luminosity, low surface-brightness, low H I mass galaxies selected 
from the ALFALFA catalog as part of The Survey for H I in Extremely Low-Mass Dwarfs (SHIELD). 
SHIELD is a comprehensive multi-wavelength project targeting 12 of the low est-H I mass galaxie s 
discovered in the ALFALFA survey which have a clear optical counterpart (jCannon et al.ll201ll ). 
The optical and H I properties of these galaxies (N. Haurberg et al. in preparation) make them 
prime candidates as XMDs. 


In Section [2] we describe the sample selection in more detail and include the derived photometric 
and H I properties. The spectroscopic observations are described in Section [3] and the reduction 
and measurement processes are outlined in Section 01 Section [5] describes how we determined 
abundances for the H II regions and the abundance results. The analysis and discussion of the 
results are in Section [H and our conclusions are in Section [8l 


2. SAMPLE SELECTION 


The ALFALFA survey includes hundreds low-H I-mass objects with Mhi < 10®Mq and additionally 
has provided the first r obust sample of gal axies at the very low-mass end of the H I mass function 
with Mhi < IO^Mq ( Martin et al. 2010 ). T welve of these galax ies, with IO^'^Mq < Mhi < 
IO^-^Mq, were selected for the SHIELD project (jCannon et al.ll201lh from a preliminary ALEALEA 
catalog. The galaxies selected for SHIELD were those with the lowest H I mass that had apparent 
optical counterparts in Sloan Digitized Sky Survey (SDSS) images. These galaxies represent the 
lowest mass potentially star-forming galaxies that have formed a significant number of stars. Thus, 
following the traditional luminosity-metallicity trend (or mass-metallicity trend), the galaxies in 
this sample should be some of the most metal-poor star-forming galaxies in the local universe. 

The SHIELD project encompasses a complex multi-wavelength data set includin g detailed H I 
gas m apping from multi-configuration Expanded Very Large Array® observations ( Cannon et al 


201 ih . broadband BVR and narrowband Ha imaging from the WIYN 3.5m telescope at Kitt Peak 
National Observatory® (KPNO; N. Haurberg et al. in preparation), 3.6 and 4.5^m imaging with 
the InfraRed Array Camera (IRAC) on the Spitzer Space Telescope^, Hubble Space Telescope® 
Imaging with the Wide Field Camera 3 using the F606W (~ R) and F814W (~ I) filters, and 
long-slit optical spectra with the Richey-Chretien Focus Spectrograph on the Mayall 4m at KPNO 
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(this work). 


Since these galaxies were selected from a blind H I survey, the sample does not suffer any bias toward 
high-surface brightness objects as is often present in optical catalogs. The only optical qualification 
for galaxies in the SHIELD sample is that there be some evidence of a stellar counterpart associated 
with the H I source in SDSS images. Low H I mass sources in the ALFLFA catalog with no apparent 
optical counterpart have been studied bv I Adams et al.l (j2013l i and are also of great interest as these 
may represent sources with extremely low surface brightness counterparts or starless dark matter 
halos. 

Thumbnail images showing the WIYN i?-band and Ha images for the 12 SHIELD galaxies are 
shown in Figure [H The H II regions identified from the displayed Ha images were targeted for 
the long-slit spectroscopic follow-up that is the subject of this paper. Figure [T] also shows the 
approximate slit locations overlaid on each galaxy. A comprehensive set of general and photometric 
quantities that describe the SHIELD galaxies is compiled in Tabled! Column 1 lists the galaxies in 
the SHIELD sample with the coordinates for each source listed in Columns 2 and 3. The distance 


calcu lated using HST optical imaging to determine the tip of the red giant branch (|McQuinn et al 


201J), is given in Column 4 and is used as the assumed distance for the derivation of all distance- 
dependent quantities. The absolute H-band magnitude, from the WIYN 3.5m (N. Haurberg et 
al. in preparation) images is listed in Column 5; the listed uncertainty includes the error in both 
the distance and photometric error. The B — V color and H-band central surface brightness not 
corrected for inclination effects are listed in Columns 6 and 7, also from N. Haurberg et al (in 

preparation). Columns 8 and 9 contain the 4.5//m and 3.6^m flux (in mJy) from the Spitzer Space 
Telescope observations, and Column 10 lists the stellar mass estimate derived from the Spitzer 
data (Cannon, Marshall, et al. in preparation). T he stellar masses presented in this table were 
calculated usi n g the method of lEskew et al.l ()2ni2l ). The H I masses given in Column 10 are from 
Cannon et al.l (j201llL but have been adjusted for the updated distances. Column 11 lists the star- 
formation rate (SFR) ca lculated from the Ha images (N. Haurberg et al. in preparation) using the 
standard prescription of iKennicuttI (jl998l i to convert from Ha luminosity to SFR. 


The optical images from the WIYN 3.5m telescope confirmed the optical counterparts of the 
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cooperative agreement by Associated Universities, Inc. 

®The WIYN Observatory is a joint facility of the University of Wisconsin-Madison, Indiana University, Yale 
University, and the National Optical Astronomy Observatory. 

^This work is based in part on observations made with the Spitzer Space Telescope, which is operated by the Jet 
Propulsion Laboratory, California Institute of Technology under a contract with NASA. Support for this work was 
provided by NASA. 

®Based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope 
Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA 
contract NAS 5-26555. These observations are associated with program 12658. 
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SHIELD H I sources as low surface brightness dirrs with blue colors (Table [T]). The Ha imag¬ 
ing revealed that the star formation in these galaxies was generally confined to 1 or 2 knots of 
star formation which were targeted with long-slit spectroscopic follow-up observations presented 
in this work. Two of the SHIELD galaxies, AGC 749241 and AGC 748778 were not included in 
this spectroscopic follow-up as the Ha imaging revealed no star-forming nebulae. AGC 749241 had 
no detected Ha emission associated with it and AGC 748778 displayed only very weak, diffuse Ha 
emission (2.6 a above the sky) with no discrete H II regions (see Figured]). Additionally, H II 
regions were detected in AGC 174585 and AGC 174605 but they were too faint to provide useful 
spectra so they are not included in the following analysis. Hence, we were able to derive useful 
abundance estimates for 8 of the 12 SHIELD galaxies. 


3. OBSERVATIONS 
3.1. OPTICAL IMAGING 

Optical imaging observations were performed using the Mini-Mosiac Imager on the WIYN 3.5m 
telescope at Kitt Peak National Observatory (KPNO) over four nights: 2010 October 7-8 and 
29-30 in March 2011. The two Fall nights (Oct. 2010) had fantastic seeing 0.5") and were 
done under photometric conditions, while the two nights in the Spring (Mar. 2011) had degraded 
seeing (~ 1.2") and were not photometric. Short-exposure post-calibration observations for the 
Spring targets were performed in April of that 2011. The broadband Johnson B, V, R, and W036 
narrowband Ha filters were used. The nominal field of view was 9.6'x 9.6' across 2 chips and 4 
amplifiers. 

Exposure times for the broadband images in the fall were 900, 720, and 600 s for B, V, and R bands, 
respectively. The narrowband imaging was taken in sets of two long narrowband Ha exposures (900 
s) with a short (180 s) R-broadband image between them. Since the conditions were less ideal in 
the Spring, longer exposure times were used. The Spring broadband data were taken with 1200, 
720, and 900 s exposures (for B, V, and R, respectively) and 1200 s narrowband Ha exposures 
with a 240 s R-band exposure in between. 

The broadband observations times should be sufficiently long to have detected any diffuse emission 
from the galaxies. These optical broadband images were analyzed with both isophotal fitting 
techniques as well as large aperture photometry as discussed in detail in N. Haurberg et al. in 
preparation. While we can not rule out that there may be an extraordinarily low surface brightness 
component that was not detected in our images, comparison of with the subsequent HST imaging 
(?????) shows now indication of such, thus we remain confident our images were sufficient at 
detecting all of the optical emission in these galaxies. 
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Fig. 1.— Thumbnail images of the SHIELD galaxies, oriented so that North points up and East 
points to the left. On the left of each column are short exposure i2-band images and on the right 
are continuum subtracted Ha images from the WIYN 3.5m (N. Haurberg et al. in preparation). 
Each thumbnail is approximately 50" x 50". The approximate position of the slit is indicated for 
the 8 SHIELD galaxies for which we obtained usable spectra. Those objects with no slit image 
superimposed either had no visible Ha emission or did not produce usable spectra. 





























































Table 1:: SHIELD Galaxy Properites 


Galaxy 

R.A. 

[2000.0] 

Dec 

[2000.0] 

Distance^ 

[Mpc] 

Ml 

[mag] 

B-V^ 

[mag] 

^0,B 

[mag "-2] 

13.63 

[mjy] 

14.53 

[mJy] 

log(M*)3 

[Mq] 

log(Mif/)"‘ 

[Mq] 

log(SFR)2 

[M0yr-l] 

AGC 748778 

00:06:34.3 

15:30:39 

6.46 ±0.15 

-10.34 ±0.08 

0.22 ±0.03 

23.78 

0.0453 

0.0310 

5.82 

6.72 

-4.92“P’ 

AGC 112521 

01:41:07.6 

27:19:24 

6.58 ±0.18 

-10.59 ±0.08 

0.45 ±0.03 

24.55 

0.563 

0.402 

6.93 

7.11 

-4.10±0.08 

AGC 110482 

01:42:17.4 

26:22:00 

7.82 ±0.21 

-13.02 ±0.13 

0.49 ±0.02 

22.45 

1.05 

0.727 

7.35 

7.31 

-2.67 ±0.04 

AGC 111946 

01:46:42.2 

26:48:05 

9.02 ±0.25 

-11.87±0.12 

0.30 ±0.02 

24.06 

0.235 

0.164 

6.83 

7.11 

-3.10±0.07 

AGC 111977 

01:55:20.2 

27:57:14 

5.96 ±0.10 

-12.60±0.09 

0.48 ±0.02 

22.67 

1.04 

0.669 

7.11 

6.85 

-3.06 ±0.08 

AGC 111164 

02:00:10.1 

28:49:52 

5.11±0.07 

-11.16±0.10 

0.42 ±0.02 

23.83 

0.470 

0.355 

6.64 

6.64 

-3.40 ±0.05 

AGC 174585 

07:36:10.3 

09:59:11 

7.89 ±0.19 

-11.32±0.13 

0.41 ±0.04 

23.17 

0.139 

0.0944 

6.48 

6.94 

-3.21 ±0.05 

AGC 174605 

07:50:21.7 

07:47:40 

10.89 ±0.28 

-12.22±0.11 

0.39 ±0.03 

22.29 

0.207 

0.147 

6.94 

7.19 

-3.19±0.04 

AGC 182595 

08:51:12.1 

27:52:48 

9.02 ±0.28 

-12.70±0.13 

0.52 ±0.03 

22.31 

0.622 

0.408 

7.25 

6.92 

-2.65 ±0.07 

AGC 731457 

10:31:55.8 

28:01:33 

11.13±0.18 

-13.73±0.11 

0.38 ±0.03 

21.31 




7.09 

-2.52 ±0.07 

AGC 749237 

12:26:23.4 

27:44:44 

11.62±0.18 

-14.12±0.12 

0.45 ±0.03 

21.53 

0.837 

0.544 

7.60 

7.49 

-2.34 ±0.05 

AGC 749241 

12:40:01.7 

26:19:19 

5.62 ±0.15 

-10.25 ±0.19 

0.14 ±0.03 

23.37 




6.74 

-5.00“P- 


^McQuinn et all l|2014h derived from HST images, tip of the red giant branch; ^N. Hanrberg et al (in pr eparation) derived from WIYN 3.5m observations; 
^ Cannon, Marshall, et al. in prep.; derived from Spitzer Space Telescope Images; ‘ ICannon et al.l ll201lh . derived from radio observations and adjusted for 
new distance estimate 














3.2. SPECTRAL OBSERVATIONS 


The spectral observations presented in this work were carried out using the Mayall 4m telescope at 
Kitt Peak National Observatory with the Richey-Chretien Focus Spectrograph and T2KA imager 
over the course of 5 nights: 2012 April 17-19 and 2012 October 15-16. The KPC-lOA grating 
(316 lines mm“^) and WG-345 blocking filter were used. The grating is blazed at 4000 A giving a 
dispersion of 2.78 A pixel”^ and total coverage from 2850-8550 A on the CCD. All spectra were 
taken with a slit-width of 1.5^^ and the slit extended 324" along the spatial direction. 

Target sources were too dim to be seen with the acquisition cameras, so nearby bright stars were 
used for blind-offsets. Offsets were checked by using multiple stars in the field, and were successful 
in aligning sources within the slit. The slit was positioned as nearly along the parallactic angle 
as possible to avoid the effects of differential refraction through the atmosphere. Since the star¬ 
forming regions in these galaxies are relatively isolated, usually only the target source fell in the 
slit, but in two cases an additional emission nebula ended up in the slit as well. 


We observed several spectrophotometric standard stars throughout each night, which were used to 
calibr ate t he flux scale for our spectra. The standards were selected from the lists of lOke fc Cunn 
(j 19831 ) and iMassev et al.l (| 19881 ). Each night we additionally took images of HeNeAr spectral lamps 
for wavelength calibration, zero-length exposure bias images, internal quartz lamps for flat field 
calibration, and twilight sky images to correct for slit-width variations. 


4. SPECTRAL REDUCTION AND MEASUREMENT 


Spectral images were processed through the standard spectral reduction routines in iraf®. The 
bias level was determined from the overscan region in each image and the mean bias image was 
used to account for any two-dimensional (2D) structure in the bias. Median combined quartz lamp 
flats were used to account for pixel-to-pixel variations and twilight sky ffats were used to create 
an illumination function correcting for variations in the slit width. The science images w ere then 
processed through the lacos_SPEC cosmic ray rejection routine of Ivan Dokkuml (j200l|) and the 
resultant “cleaned” images were examined carefully by eye in comparison with the original image 
to ensure that emission line pixels were not rejected as cosmic rays. 


H II region spectra were extracted using the apall package in iraf. Since the diffuse background 
emission was typically very weak, local sky regions that closely bracketed the emission line source 
were selected for night-sky subtraction. The wavelength scale was determined for each night using 
the solution derived from HeNeAr lamp spectra. The spectrophotometric standard stars were then 
used to create a sensitivity function that was applied to all the spectra. After the spectra were 


®IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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wavelength and flux calibrated, standard stars were compared to a catalog spectrum in order to 
check the calibration, which appeared good in all cases. As a last step, telluric absorption features 
were removed. 


Figure [2] shows examples of the final spectra. Several important strong lines were present in all 
the extracted spectra, most notably Ha, H/3, the [O II] doublet, and the strong lines from the 
[O III] triplet. The weak, temperature sensitive, [O III] A4363 line was discernible in two of the H II 
regions, AGC 1104826 and AGC 182595. The other extracted spectra were good quality but did 
not contain an unambiguous detection of the [O III] A4363 line. In total, we extracted 11 unique 
H II regions from 8 galaxies. Both of the H II regions with discernible [O III] A4363 are included in 
Figure [2j We have plotted a version showing the full dynamic range of the emission and another 
version that displays only a small range in flux. 

Once emission lines were identified they were measured using the splot routine in iraf. This 
routine allows interactive setting of the continuum level at each line. To decide on the correct 
location for the continuum we used an average of the local continuum on either side of the emission 
line and assumed a linear fit for for the continuum under the emission line. We generally assumed 
that the continuum remained at a constant level under the emission line except for those cases where 
the lines were very broad or the slope of the continuum was steep enough to change significantly 
through the line. 


Since we are measuring emission from regions of relatively recent star-formation it is expected that 
the permitted Balmer lines are affected by underlying absorption from hot main-sequence stars. The 
underlying absorption was accounted for by first assuming that the equivalent width of absorption 
is the same for all the Balmer lines and then using multiple Balmer emission-line ratios (Ha/H/3, 
H7/H/3, and H(5/H/3) to simultaneously estimate the amount of underlying absorption and c(H/3) 
(i.e. the absorption coefficient at H/3). We were able to estimate the proper absorption correction 
by determining where the c(H/3) coefficients calculated from each Balmer ratio converged. This 
correction was applied to the hydrogen lines so that the relative Balmer intensities could be used 
to accurately determine the reddening for each H II region. In cases where only one Balmer ratio 
was available, an average absorption correction of 1.0 A was used since most H II regions showed 
some absorption th at we felt need e d to be accounted for. The reddening was determined using 
the reddening law of iGardelli et al.l (119891 ) and using the temperature derived from the [O III] lines 
when possible. If the [O III] A4363 line was not detected, we assumed an electron temperature of 
10,000 K. 


Since the Balmer decrement was used to calculate the reddening, there is no distinction between 
intrinsic and foreground reddening. The reddening corrected flux measurements (relative to H/3) 
for each line we measured are tabulated for each H II region in Table O along with the as¬ 
sumed equivalent width of underlying absorption (Wabs), c(H/3), and the assumed Tg and ng. 
The electron temperature (Tg) and electron density (ne) of the nebula were estimated using the 
[O III] A4363/[0 III] AA4959,5007 and [S II] AA6716,6731 ratios, respectively. For cases where the 
[O III] A4363 line was measured, the temperature and density calculation were performed with the 
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ELSA program (| Johnson et al.ll2006l l described more completely in the following section. In cases 
where there was not a reliable measnrement of [S II] A A6716,6731 a default electron density of 100 
cm“^ was assumed. This density is consistent with what is commonly observed for H II regions. 
The errors given on the line fluxes were calculated by fully propagating the derived errors from the 
relevant quantities (e.g., rms of the continuum, rms scatter in flux calibration, etc.). 


5. ABUNDANCE DETERMINATION 


The emission line ratio of [N II] A6583/Ha and [O III] A5OO7/H;0 can be used to plot each H I l-region 


on a traditional diagnostic diagram (e.g., [Baldwin et al. 


1981 


Veillenx fc Osterbrocklll987l ) shown 


in Figure [3l which provides a rough diagnostic for excitation and temperature (thus metallicity) 
of emission nebulae. SHIELD galaxies are plotted on Figure [3] as solid black points; those with 
circles around them are the two H II regions which displayed a measurable [O III] A4363 flux. The 
green triangles represent H II regions in low H I mass dirr galaxies from t he Arecibo Dual Beam 
Survey (ADBS) which is a similarly selected, but more luminous, sample (jHaurberg et al.ll2013l L 
and the grey poin t s are galaxies from the KPNO International S pectroscopic Survey (KISS) (e.g.. 


Salzer et al 


200C; 


Melbourne fc Salzer 


2002 ; 


Salzer et al 


20051) . The dashed line differentiates 


active galactic nuclei (AGN) from star-forming nebulae (iKanffmann et al.ll2003l ) and the solid line 
repres ents the locus of high excitation star-forming nebulae from the models of iDopita fc Evans 
(|l986l ). which increase smoothly in metallicity from the upper left to the lower right. 


It is clear that H II regions from the SHIELD galaxies lie below and to the left of the majority of 
the KISS sample and the theoretical high-excitation star-forming model curve. This is consistent 
with the manner in which these galaxies were selected and indicates that the observed H II regions 
in the SHIELD galaxies are generally low-excitation systems. Emission-line surveys like KISS are 
much less likely to include these types of low-excitation systems because the key emission lines used 
for selection in those surveys need to be strong, high equivalent-width lines in order to be detected. 
The relatively low-excitation of these systems indicates that either the initial mass function (IMF) 
was not fully populated or the most massive stars formed in these regions have already evolved off 
the main sequence. Both possibilities are consistent with the expectation for very low-luminosity 
galaxies as either would lead to reduced luminosity and surface brightness in comparison to the 
more well-studied low-metallicity BCDs. 


Three of the H II regions that we observed (AGO II04826, AGG II252I, AGC III977) had such 
weak [N II] A6583 lines that we were unable to determine an accurate measurement of [N II] A6583/Ha 
Thus, these three H II regions have been plotted with arrows indicating the measured [N II] A6583 /Hq; 
is an upper limit. The lack of distinguishable [N II] A6583 emission is consistent with these being 
among the lowest metallicity regions (see Table [3]). The H II regions from the SHIELD galaxies all 
occupy the region of the diagram expected for low-metallicity star-forming systems. However, the 
H II regions spread over a significant portion of the low-metallicity region while th e derived abun 


dances in Table [3] cover a quite narrow range. This is consistent with the finding of lHaurberg et al 
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Fig. 2.— Four examples of spectra from H II regions in the SHIELD galaxies. The first two display 
a measurable [O III]A4363 emission line and the location of this line is marked with a small arrow. 
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Fig. 3.— Emission-line diagnostic diagram with each H II region plotted separately. The black 
points are from the SHIELD galaxies in this work, those with outer circles are H II regions for which 
we obtained a measurement of th e temperature sensitiv e [O III] A4363 line. The green triangles are 
galaxies from the ADDS survey ([Haurberg et al.l 1201311 and the grey points are a large sample o f 
emission line sources from the KISS survey (e.g., ISalzer et al.l l2000l: [Melbourne &: Salzed l2002l L 


The solid line traces th e star-forming galaxy models of 


Dopita fc EvansI "i 1986): 


the dashed line 


(jKauffmann et al.l 120031 ) delineates the region of the diagram usually occupied by AGN. 



























Table 2;: Line Measurement Results 


Line Identification 

AGO 112521 

AGO 110482 a 

AGC 1104826 

AGC 111946 

AGCl 11977 

AGC 111164 

(A) 

Fion/H/I 

Fion/H/? 

Fion/H/3 

Fion/H/3 

Fion/H/3 

Fion/H/3 


[0 II] A3728 

0.809 

± 

0.075 

3.062 

± 

0.217 

1.306 

± 

0.091 

3.137 

± 

0.282 

1.965 

± 

0.250 

1.417 

± 

0.095 

[Ne III] A3869 




0.300 

± 

0.031 

0.555 

± 

0.044 




1.161 

± 

0.150 

0.213 

± 

0.019 

He I+HC 
















0.152 

± 

0.015 

[Ne III] A3970+He 







0.077 

± 

0.023 







0.133 

± 

0.013 

m 

0.260 

± 

0.027 

0.432 

± 

0.034 

0.252 

± 

0.027 

0.298 

± 

0.045 




0.249 

± 

0.018 

H7 

0.387 

± 

0.032 

0.433 

± 

0.031 

0.520 

± 

0.034 

0.508 

± 

0.052 




0.436 

± 

0.025 

[0 III] A4363 







0.109 

± 

0.021 










H/3 

1.000 

± 

0.060 

1.000 

± 

0.051 

1.000 

± 

0.045 

1.000 

± 

0.068 

1.000 

± 

0.164 

1.000 

± 

0.044 

[0 III] A4959 

0.527 

± 

0.036 

0.939 

± 

0.049 

1.640 

± 

0.066 

0.276 

± 

0.034 

0.857 

± 

0.090 

0.750 

± 

0.035 

[0 III] A5007 

1.701 

± 

0.092 

2.708 

± 

0.124 

5.189 

± 

0.189 

0.838 


0.061 

3.375 

± 

0.270 

2.339 

± 

0.096 

He IA5876 

0.113 

± 

0.017 




0.051 

± 

0.014 










Ha 

2.860 

± 

0.225 

2.860 

± 

0.201 

2.785 

± 

0.169 

2.860 

± 

0.237 

2.860 

± 

0.510 

2.860 

± 

0.183 

[N II] A6583 

0.015 

± 

0.012 

0.080 

± 

0.014 




0.226 

± 

0.029 

0.176 

± 

0.044 

0.062 

± 

0.008 

[S II] A6717 

0.143 

± 

0.020 

0.249 

± 

0.024 

0.157 

± 

0.016 

0.227 


0.029 




0.128 

± 

0.012 

[S II] A6731 

0.131 

± 

0.019 

0.179 

± 

0.019 

0.124 

± 

0.015 

0.153 

± 

0.024 




0.094 

± 

0.010 

[Ar III] A7136 







0.060 

± 

0.012 










Te 

10,000 

K* 

10,000 

K* 

15, 

63C 

K 

10,000 

K* 

10,000 

K* 

10,000 

K* 

He 


389 



48 



146 

100** 

100** 


71 


Wabs (A) 


5.0 



2.5 



1.0 



0.0 



1.0 



2.5 


c(H/3) 

0 

.069 

0 

.277 

0 

.496 

0 

.532 

0 

.375 

0 

.171 


Continued on next page 
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Line Identification 

AGO 182595 

AGC 731457 

AGC 749237a 

AGC 7492376 

AGC 749237c 

(A) 

Fion/H/3 

Fion/H/3 

Fion/H/3 

Fion/H/3 

Fion/H/3 

[0 II] A3728 

1.941 

± 

0.183 

3.705 

± 

0.353 

3.283 ± 

0.365 

3.326 

± 

0.314 

3.212 ± 0.354 

[Ne III] A3869 

0.493 

± 

0.048 






0.622 

± 

0.060 

0.475 ± 0.062 

He I+HC 













[Ne III] A3970+He 













m 

0.256 

± 

0.026 

0.138 

± 

0.028 



0.226 

± 

0.026 


Hy 

0.429 

± 

0.036 

0.319 

± 

0.036 

0.421 ± 

0.051 

0.435 

± 

0.038 


[0 III] A4363 

0.087 

± 

0.015 










H/3 

1.000 

± 

0.067 

1.000 

± 

0.069 

1.000 ± 

0.085 

1.000 

± 

0.070 

1.000 ± 0.086 

[0 III] A4959 

1.373 

± 

0.091 

0.593 

± 

0.047 

0.892 ± 

0.079 

0.796 

± 

0.057 

0.709 ± 0.065 

[0 III] A5007 

4.159 

± 

0.264 

1.784 

± 

0.120 

2.487 ± 

0.194 

2.102 

± 

0.139 

2.118 ± 0.165 

He IA5876 













Ha 

2.786 

± 

0.264 

2.860 

± 

0.270 

2.860 ± 

0.315 

2.860 

± 

0.277 

2.860 ± 0.319 

[N II] A6583 

0.064 

± 

0.013 

0.143 

± 

0.023 

0.197 ± 

0.036 

0.128 

± 

0.020 

0.120 ± 0.029 

[S II] A6717 

0.252 

± 

0.029 

0.210 

± 

0.029 

0.186 ± 

0.036 

0.446 

± 

0.050 

0.313 ± 0.047 

[S II] A6731 

0.193 

± 

0.024 

0.160 

± 

0.025 

0.179 ± 

0.035 

0.319 

± 

0.037 

0.132 ± 0.030 

[Ar III] A7136 

0.085 

± 

0.015 










Tg 

15,680 K 

10,000 K* 

10,000 K* 

10,000 K* 

10,000 K* 

ng 


127 


119 


489 



41 


100** 

Wabs (A) 


2.0 



5.0 


0.0 



1.0 


1.0“^ 

c(H/3) 

0.191 

0.286 

0.044 

0.051 

0.074 


** indicates a default Ug of 100 cm“^ was assumed. Tg values marked with * were not 
calculated directly, but were assumed to be 10,000 K. next to the value for Wabs 
indicates that the Balmer line data was not sufficient to derive a value for Wabs so an 
average absorption correction of 1.0 A was used instead. 
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(|2013l l that there is significant scatter in the position of similar metallicity H II regions on the 
diagnostic diagram. 

In order to accurately determine the chemical abundances for emission line nebulae, it is preferable 
to have a direct measurement of the electron temperature Tg. In the wavelength range covered 
by our observations, the best way to calculate Tg is from measuring the [O III] A4363 line. The 
[O III] A4363 and [O III] A A4959, 5007 lines are produced by electronic transitions from different 
energy levels in the ion. The relative strengths of the lines thus depend on the population of 
the different energy levels, which is strongly dependent on the temperature of the electrons (and to 
a lesser extent, the electron density). [O III] A4363 is usually a very weak line that is often difficult 
to detect. We were able to detect this line in two of our H II regions, although the detection in AGC 
182595 is noisy and thus fairly uncertain. In order to calculate the nebular abundances for these 
H II regi ons we used th e ELSA (Emission Line Spectrum Analyzer) program which is described in 
detail in l.Iohnson et all (2006). 

ELSA uses a five-level atom routine with ionization correction factors and a two-region ioniza- 
tion model t o calc ulate abundances. The five-level atom calculations are based on the w ork of 
Henry et al. ( lOsd l but include multiple updates and improvements (see Johnson et al. 2006). The 


[O III] A4363/[0 III] A A4959,5007 ratios is used to calculate temperature in the high-ionization re¬ 
gion and the [S II] A A6716, 6731 ratio used to calculate the e lectron densi t y (in the low-ionization re¬ 
gion) using an iterative process similar to that described in Izotov et al. ( 2006 I l The low-ionization 
region temperature would ideally be calculated from the ratio [N II]A5755/[N II]AA6548,6583, how¬ 
ever, [N II]A5755 is a weak line that is usually not detectable in star-forming regions. Therefore, 
the temperature o f the low-ionization region is estimated using Tg q +2 and following the method of 
Pagel et al. 1 19921 ). The ionization correction factors are determined from lines with similar ion¬ 
ization potentials and then used to account for the unseen ionization states. Many uncertainties, 
including those in line fluxes, reddening correction, plasma diagnostics, and ionic abundances are 
carried through the ELSA program and propagated properly into the final uncertainties in the 
abundances. 


In cases where Tg is not directly measurable, strong-line calibrations must be used. These methods 
rely on ratios of strong emission lines to estimate an electron temperature and thus metallicity. 
Specifically, most strong-line abundance methods are calibrated to a specific oxygen abundance; 
we will use oxygen abundance to be synonymous with metallicity for the remainder of this paper. 
While these methods are very useful, some caution must be employed as they are reliant on both 
theoretical models and empirical calibration to arrive at an abundance. 

One of the most robust strong-line methods for determining oxygen abundance is that of McGaugh 
(1991; henceforth the McGaugh method). The McGaugh method uses several key abundance 
ratios: 


R 23 : log(([OlIl]AA4959,5007+[O Il]]A3727)/H/3) 
O 23 ; log([0 III] A A4959,5007/[O II] A3727) 
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[N II] / [O II] : log( [N II] A6583/ [O II] A3727). 

Each measured H II region is placed on a grid of model emission nebula using the R 23 and O 23 
parameters as shown in Figure |4| Since the models are degenerate (i.e. there are two possible 
oxygen abundances that can produce the same values for R 23 and O 23 ) the value of [N 11] /[O II] is 
used to break the degeneracy between high- and low- metallicity “branches.” H II regions are placed 
on the low-metallicity branch if [N n]/[0 II] < —1.0 and the high-metallicity branch if [N n]/[0 II] 
> —0.9. If —1.0 < [N Il]/[0 II] < —0.9 the H II regions are deemed to be in the “turn-around 
region” where precise determination of the abundance can be tricky. However, all the H II regions 
in this work had very low [N II] A6583/[0 II] A3727 placing them securely onto the low-metallicity 
branch. 


The degeneracy of the McGaugh models can be troublesome and lead to significant uncertainty 
when the branch is not well determined. However, the method is ideal for use with a sample like 
ours where the branch is unambiguous. The McGaugh method is a particularly robust strong-line 
calibration because it uses the total flux in the strong oxygen lines (the R 23 parameter) to determine 
oxygen abundance, and also accounts for different hardness of the ionizing spectra by incorporating 
the O 23 parameter which greatly improves its accuracy. The precise value of the oxygen abundance 
derived from the McGaugh method depends not only on the values of R 23 and O 23 but also on 
the specihcs of w hich models were used to create the grids. We have used the models described in 


McGaugh! (|l99lh where the highest mass star produced from the IMF is 60 Mq. The exact value of 


the oxygen abundance was determined by using the position of each H II region and interpolating 
between the model points to estimate abundance more precisely. 


A final table of abundances as well as key line ratios is presented in Table |3j Golumn 1 lists the 
galaxy and specific H Il-region where applicable. Golumns 2—4 contain the line ratios used to cal¬ 
culate the abundance via the McGaugh method. Golumn 5 gives the oxygen abundance determined 
via the McGaugh method (in terms of log(0/H)-|-12) while Golumn 6 gives the oxygen abundance 
determined via the Tg method when available. Column 7 lists the error-weighted mean oxygen 
abundance from the individual H II regions in AGC 110482 and AGC 749237. The abundance 
listed in bold is the assumed abundance for each galaxy. We do not give explicit errors on the 
McGaugh abundances, but will assume an error of 0.10 dex. 


The abundances derived via the McGaugh and Tg methods are not in perfect agreement with each 
other which is consistent with what has been reported by other authors. The inconsistency is often 
considered an effect of intrinsic errors in the strong-line calibration methods since they rely on 
indirect determinations of Tg. While the Tg method provides a direct measurement of the electron 
temperature it is also subject to some level of uncertainty and bias as it has been shown to favor 
higher ten iperature region s in the case that there are temperature fluctuations in the nebula being 
observed (jPeimbertI Il967l ). Thus, it can be difficult to determine which method more accurately 
yields the “correct” abundance. 

Some authors (e.g., Kewlev h Ellisonll2008 : Zahid et al.ll2012l) have suggested that empirical calibra- 
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Fig. 4.— Theoretical grids for estimating abundance from iMcGaughl (jl99ll ) plotted with H II 
regions from the SHIELD galaxies. The roughly horizontal lines indicate different ionization pa¬ 
rameters and the roughly vertical lines represent lines of constant oxygen abundance (expressed as 
log(0/H)-|-12). A line of constant oxygen abundance is displayed for each 0.1 dex increase from 
log(0/H)-|-12 = 7.4 to 9.3. All the objects in this study appear to be located on the “far-side,” 
or low-metallicity branch (log(0/H)-|-12 < 8.2), of the surface traced by the models. The points 
plotted within circles represent those H II regions where the [O III] A4363 line was detected. 

























Table 3:: Derived SHIELD Galaxy Abundances 


Galaxy 

l0g(R23) 

log(023) 

log([NH]/[OH]) 

log(0/H)+12 

]McGaugh] 

log(0/H)-bl2 

[Te] 

log(0/H)±12 
[w. mean] 

ACC 112521 

ACC 110482 

0.482 

0.440 

< -1.732 

7.41 


7.81 ±0.09 

a 

0.827 

0.076 

-1.583 

8.05 


b 

0.910 

0.718 

< -1.685 

7.98 

7.74 ±0.08 


ACC 111946 

0.629 

-0.450 

-1.142 

7.94 



ACC 111977 

0.792 

0.333 

< -1.048 

7.88 



ACC 111164 

0.654 

0.338 

-1.359 

7.67 



ACC 182595 

0.874 

0.455 

-1.482 

7.99 

7.75 ±0.13 


ACC 731457 

ACC 749237 

0.784 

-0.193 

-1.413 

8.08 


8.03 ±0.09 

a 

0.817 

-5.30e-4 

-1.222 

8.06 


b 

0.794 

-0.060 

-1.415 

8.03 



c 

0.781 

-0.055 

-1.428 

8.01 




Calculated abundances for SHIELD galaxies. When more than one HII region was present an error weighted 
mean, listed in Column 7, was calculated. The abundances listed in bold are the abundances we assumed for 
each galaxy. The < indicates the [Nil] line was not measurable so the [Nil]/[OH] in this table is only an upper 
limit. 
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tions should be made to correct for systematic differences between various abundance determination 
methods. With only two Tg abundances, we are unable to gauge whether we see a systematic offset 
in our data. However, i f we u se the sample of low-luminosity Local Volume Legacy (LVL) galaxies 
analyzed in lBerg et al.l (|2012l L which we will use as a comparison set in Section [6l we can determine 
a reasonable empirical correction a nd use it to reconc ile the two abundance methods. We have as¬ 
sumed the Tg abundances listed in iBerg et al.l (120121 ) then calculated McGaugh abundances from 
the fluxes given in that paper, using the same method and models used for our galaxies. When the 
two different abundance calculation methods are compared we can see that there does appear to 
be a systematic offset. 


A plot of abundance diffe rence (calcu l ated a s the McGaugh abundance minus the Tg abundance) 
vs. Tg abundance for the iBerg et al.l (120121) sample is s hown in Figure [5j The red diamonds are 
H 11 regions from the the LVL galaxies in iBerg et al.l ()2012l i and the black points are the H 11 
regions with Tg abundances from the SHIELD galaxies. Since it is unclear whether this difference 
is a constant offset or if it varies with metallicity we have only used H II regions within the range 
of metallicities that is consistent wit h our samp l e (7.4 < log(0/H)-|-12 < 8.1) and exclude the 
obvious outlier. The points from the iBerg et al.l (120121 1 sample which were used to calculate the 
mean offset are plotted as filled red diamonds and those that were excluded are plotted as open 
red diamonds. The mean offset that we calculated was 0.076 dex, indicating that on average the 
McGaugh abundances were 0.076 dex higher than those calculated using the Tg method This offset 
will later be used as a “correction” when comparing our sample to others from the literature. 


“Corrected” average abundances for the SHIELD galaxies are compiled in Table IH When more than 
one H II region was measured a weighted mean of the multiple H II region abundances was used. 
Notice that the abundances for AGC 112521 and AGO 111164 qualify them both as XMD galaxies 
which is consistent with them being the two lowest luminosity galaxies for which we obtained usable 
spectra. 


6. ANALYSIS 

There is a well-established relationship between the luminosity and metallicity of galaxies, where 


luminosity dirrs (e.g.. 

Leoueux et al.lll97f 

): Skillman et al. 

1989a 

1: PilvuginIboOl: 

Melbourne & Salzer 

2002; 

Tremonti et al. 

2004; 

Salzer et al.l 

>nn5bl: 

Lee et al. 

2006; 

van Zee & Havnes 2006: 

Berg et al. 

2012; 

Haurberg et al. 

2013 

Skillman et al. 

2013'). The underlving cause of this relationship is not 


fully understood, but it is believed to trace a more fundamental mass-metallicity (M-Z) relation¬ 
ship, where luminosity is representative of stellar mass. Studies of the luminosity-metallicity (L-Z) 
relationship that use near-in frared (NIR) luininosities (which is better correlate d with stellar mass) 
show decreased scatter fe.g.. lLee et al.ll2nr)4l : ISalzer et ahlDDOhbl : iLee et al.lPOObl ). implying the M-Z 
relationship is the more fundamental physical relationship. The origin of the relationship, however, 
is still unclear but is a key to understanding the star-formation histories and chemical evolution of 
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Log(0/H) + 12 [Tj 


Fig. 5.— Comparison of derived McGaugh and Tg abundances. Black points within circles are H II 
regions from galaxies in t h e SHI ELD sample with Tg derived abundances. Red diamonds represent 
galaxies from iBerg et al.l (j2ni2l ) ; filled diamonds are those used to calibrate the empirical offset 
between the two abundance methods, open diamonds are either obvious outliers or outside of the 
abundance range that we considered for the offset calculation. 











































Table 4:; “Corrected” SHIELD Galaxy Abundances 


Galaxy 

log(0/H)+12 

Note 

ACC 112521 

7.33 ±0.10 

corrected McGaugh, XMD galaxy 

ACC 110482 

7.79 ±0.07 

weighted mean of a (corrected) and b (Tg) 

ACC 111946 

7.86 ±0.10 

corrected McGaugh, 

ACC 111977 

7.80 ±0.10 

corrected McGaugh 

ACC 111164 

7.59 ±0.10 

corrected McGaugh, XMD galaxy 

ACC 182595 

7.75 ±0.13 

Tg abundance 

ACC 731457 

8.00 ±0.10 

corrected McGaugh 

ACC 749237 

7.95 ±0.06 

weighted mean of a, b, and c (all corrected) 


Corrected abundances for SHIELD galaxies (see text for correction). For ACC 110482 and ACC 
749237 a weighted mean of the abundances for each HII regions was used. Both ACC 112521 and 
ACC 111164 qualify as XMD galaxies. 
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dwarf galaxies. 


Low luminosity galaxies (e.g., Mb > —13), which should also be the most metal-poor galaxies, 
are a relatively understudied population, thus the nature of the L-Z and M-Z relationships in the 
very low-metallicity regime remains indeterminate. The XMD galaxies which have been studied are 
mostly BCDs with very high star-formation densities and relatively high star-formation rates and 
are not believed to be representative of the majority of the very low-mass dirr population. A much 
more in-depth analysis of extremely low-metallicity systems, including low surface brightness dirrs, 
is necessary to characterize these relationships. In Figure [6] we have plotted the SHIELD sample 
in luminosity-metallicity diagrams using oxygen abundance (plotted as log(0/H)-|-12) to represent 
metallicity with Mb on the luminosity axis. On Figure [Gjo and c we have also included comparison 
galaxies from the literature, described in more detail in the paragraph below. The SHIELD galaxies 
(black points) within an outer circle represent galaxies where a Tg abundance was measured. The 
Te abundance in Table [3] has been adopted for AGO 1104826 and the weighted mean of the 3 H II 
regions in AGC 749237 has been adopted for that system. 


Our comparison sample is described in Table [5j It has been broken into two pieces: the “primary 
sample (plotted as filled symbols) i s compose d of o nl y those galaxies w it h vel ocity-i ndependent 


distan ces and Tg abundances from iLee et al 


( 2006l i: Ivan Zee fc HavnesI ( 20061). and Berg et al 


(I 2 OI 2 IL The “expanded ” samp le (plotted as open symbols) includes galaxies from iBerg et al.l (|2012l i 
and Ivan Zee fc HavnesI (1200611 with velocity-b ased distances or strong-line abundances as well as 
the ADBS sample from iHaurberg et al.l (120131 ). The ADBS sample is of particular interest as it was 
chosen in a similar manner to the SHIELD sample (e.g., ADBS is a blind H I survey) and thus has 
similar selection biases. The ADBS galaxies have velocity-derived distances and the abundances 
were derived using the same techniques presented in this work with a mix of both the Tg and 
McGaugh methods used to derive the abundances (see Haurberg et al. 2013 for details). 


Figure [6^ shows just the SHIELD galaxies, while Figure [6)3 and c include comparison samples from 
the literature. Figure [6)3 includes the primary and expanded comparison samples while Figure [6]c 
shows only the primary comparison sample and the SHIELD galaxies. As indicated, the abundances 
from the SHIELD galaxies in Figure [6b have been corrected, meaning the McGaugh abundances 
have been adjusted by the offset calculated in Section [5] in order to be consistent with the Te 
scale. 


The data shown in Figure [6] were fit using a bivariate linear regress ion method weighting the e rrors 
in both X and y. The method we used is based on that outlined in lAkritas &: Bershadvl (Il996l l and 
references therein. The fit to the primary comparison sample is shown on all three panels in Figure 
[6] as a solid line. The dash-dot line shown in Figure [6] b is the fit to all of the comparison sample 
plus the SHIELD galaxies and the dashed line on Figure [6b is the fit to the primary comparison 
sample plus the SHIELD galaxies (using the corrected abundances). The XMD outlier, UGC 5340 
{Mb = —15.83, log(0/H)-|-12 = 7.20), was excluded from all of the fits. The fit that is shown with 
a dotted line on Figure [6)3 was obtained by fitting only the H I-selected galaxies (SHIELD and 
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Table 5:: Comparison Sample 


Sample Reference 


Symbol 


Primary 


Expanded 


Stellar based distances and Tg abundances 
Berg et al. (2012) Filled Red Diamonds 

Lee et al. (2006 1 Filled Blue Points 

van Zee Havnes (2006) Filled Orange Squares 

Velocity determined dist. and/or strong-line abundance 
Haurberg et al. (2013 1 Open Green Triangles 
Berg et al. (20121 Open Red Diamonds 

van Zee fc Havnes (2006 1 Open Orange Squares 
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ADBS samples) and is interesting since it is notably different from the other fits. The fit parameters 
and samples which were used to obtain them are described in Table [6l The reduced values listed 
in that table are calculated for each fit and compared to the sample plotted in the indicated figure; 
Xred calculated assuming abundance to be the dependent variable. We chose not to include 
a fit to the SHIELD sample alone because it is a statistically small sample that covers a limited 
luminosity range, thus it produced significant uncertainties in the quality of the fit. 

The various fits are in general agreement within the errors (with the exception of the “H I - sample” 
fit which is discussed in the following paragraph). In all cases, the inclusion of the SHIELD sample 
results in a more shallow slope, as could easily be predicted from inspection of Figure [H The 
SHIELD galaxies clearly lie above literature galaxies of similar luminosity even when the corrected 
abundances are used. This is certainly an intriguing trend, but since we have a relatively small 
sample it is possible that this is due to intrinsic scatter in the relationship owing to differing mass- 
to-light ratios in the H-band or uncertainty in measurement of both luminosity and abundance. 
However, it may represent something more significant such as a fundamental difference in the 
samples or a flattening of the L-Z relationship at low-luminosity. 

The fit obtained from the H I-selected samples (dotted line) is interesting as it is significantly 
different from the other fits. It has a slope that is substantially shallower than the average slope 
of the other fits, but provides the best fit to the SHIELD sample according the reduced test. 
This indicates that the sample selection process may be a contributing factor and may suggest 
underlying physical differences between samples selected from H I surveys and those selected from 
optical catalogs or emission line-surveys. For consistency, we used the uncorrected abundances 
for Figure [6)b and when fitting this H I selected sample since the ADBS galaxies have a mix of 
strong-line and Tg derived abundances. When the corrected abundances are used (Fig. [Uc) the 
obvious discrepancy between the SHIELD sample and the fit from the comparison sample is reduced. 
However, even when the corrected abundances are used, all but one of the SHIELD galaxies lie 
above the fit line indicating some level of systematic difference that is unlikely to be due solely to 
intrinsic scatter in the relationship. These issues are further discussed in Section [71 

Figure [7^ shows another version of the L-Z diagram. This version uses the corrected abundances 
for the SHIELD galaxies, like FigureEjc, but uses the NIR luminosity (4.5/im band; M 4 . 5 ) instead of 
the H-band luminosity as the luminosity metric. Since we do not have reliable 4.5^m luminosities 
for one of the SHIELD galaxies for which we measured an abundance (AGC 731457), this version of 
the L-Z relationship only features 7 of the galaxies in the SHIELD sample. The SHIELD galaxies 
are more coincident with the comparison data in the NIR version of the L-Z diagram than they are 
in Figure [ 6 l The fits derived from fitting the comparison sample (solid line) are nearly identical to 
those calculated when the SHIELD data are added to the set (dashed line). Despite the similarity 
of the fits, the majority of the SHIELD galaxies still appear above the fit line on this diagram, 
again suggesting the possible presence of a systematic difference from the comparison sample. The 
fit parameters for Figure [7| are given in Table El 

In Figure |7)o the mass metallicity (M-Z) relationship is shown. For the comparison sample, stellar 
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Fig. 6.— Luminosity-metallicity diagrams. The SHIELD sample is displayed with filled black 
points; dwarf irregular galaxies from the literature are shown with various symbols, filled points 
represent the “prim ary” comparison sample and ope n points the “expande d” comparison sam- 
ple (see Table [5]): iLee et al.l (j2006l i = blue po ints. iHaurberg et ahl (j2013l i = green triangles, 
van Zee &: HavnesI (j2006l i = orange squares, and iBerg et ahl (|2012l f = red diamonds. Figure [6)b 


shows the entire comparison sample and the uncorrected SHIELD abundances (Table [3|); c includes 
galaxies from only the primary comparison sample and the “corrected” SHIELD abundances (Table 
H]) . The lines represent fits described in the text and Table El 


























































Table 6:: L-Z and M-Z Relationship Fit Parameters 


Sample 

Slope 

Intercept 

Symbol 


Goodness of Fit (x^ed) 

log(0/H)-t-12 vs. Mb . 




Ell 

Eb 

Et 

Primary Comparison 

-0.125 ±0.008 

6.01 ±0.12 

Solid 

8.54 

13.23 

12.42 

SHIELD+ Primary + Expanded 

-0.110 ±0.006 

6.24±0.10 

Dot-Dash 


11.43 


Hl-Selected (SHIELD + ADBS) 

-0.061 ±0.012 

7.00 ±0.24 

Dotted 


12.91 


SHIELDt + Primary 

-0.113 ±0.007 

6.21±0.11 

Dashed 



11.52 







Ell 







Primary Comparison 

-0.112±0.007 

5.94±0.13 

Solid 



10.07 

SHIELDt + Primary 

-0.110 ±0.008 

5.98±0.14 

Dashed 



9.76 

log(0/H)+12 vs. log(M*) - 






m> 

Primary Comparison 

0.29 ±0.09 

5.64±0.71 

Solid 



11.42 

SHIELDt + Primary 

0.30 ±0.25 

5.60 ±1.86 

Dashed 



11.57 


The results of various fits to the L-Z and M-Z relationship. ^ indicates that corrected McGaugh abundances were 
used for the SHIELD galaxies. The “goodness of fit” parameter is the reduced value for the given fit when 
compared with the data plotted on the indicated panel or figure. See text for details of the fitting process and 
different samples. 
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mass was deriv ed using the 4.5^m flux and the m 4.5 — K (ii'-band magnitude) color following 


Lee et alJ (j2006l '). The SHIELD IR observations do not include the K-band, so stellar masses for the 


SHIELD sample had to be derived in a diffe rent manner. We used the 3.6 and 4.5/rm magnitudes and 
followed the method of lEskew et al.l (j2012l L Since there can be significant uncertainties associated 
with stellar mass calculation s, we performed an independent derivation to check against the mass 
obtained with the method of lEskew et al.l (j2ni2l ). As a check we use d a straight-forward conversion 
of 3.6/rm luminosity in to mass a,ssumi ng a mass-to-light ratio of 0.5 (jMcGaugh fc Schombertll2ni3ll 


and = 3.24 (|Oh et al.ll2008l L The difference between the mass derived using lEskew et al 


(|2012l l and that derived from our “straight-forward” conversion of 3.6/im luminosity was very small; 
the average difference w as only 0.008 dex . While we can not confirm t hat o ur mass estimates are 
fully consistent with the iLee et al.l (120061') method u sed b v iBerg et al.l (j2012lL we do feel confident 


that the stellar-mass estimates from lEskew et al.l (j2012l l for the SHIELD galaxies given in this 
paper are consistent with other IR-based mass derivation methods. Due to the use of different 
stellar-mass calculation methods and the inherently difficult nature of determining stellar masses, 
we proceed with some caution concerning the results derived using these mass estimates. 


Similar to the NIR-LZ relationship, the SHIELD galaxies appear consistent with the literature 
samples in the M-Z relationship shown in Figure [7)3 and the inclusion of the SHIELD data has 
very little effect on the results of the fit. The fit parameters and x^ed M-Z relationship are 

included in Table [H 


DISCUSSION 


Various physical mechanisms have been suggested to explain the origin of the M-Z relationship for 
dirrs, but the ejection of enriched gas fr om low-mass galaxies via supernova and stellar winds is 


one of the most commonly invoked (e.g., iTremonti et al.ll2nn4l ). In this scenario, the relationship 


arises as an effect of an enriched-gas-retention sequence where more massive galaxies are better 
able to retain enriched ejecta thus producing higher gas phase metallicities. While the retention 
of metal-enriched gas is likely a contributing factor, most observational studies that indicate the 


Papaderos et al. 

1994; 

Marlowe et al. 

1995: 

Martin 

2005; 

Oev et al. 

2007) 

and may not apply to 

dirrs with much less active modes of star-formation. 

Stinson et al. 

(200 

7), among others, have 


proposed models for the evolution of dwarf galaxies where star-formation occurs in patterns of 
relatively strong starbursts that are quenched by supernova feedback and then followed by periods 
of more moderate star formation (the so called “breathing” model). This “bursty” star forma¬ 
tion history may supply a mechanism to remove enriched gas while allowing for periods of more 
quiescent evolution. However, it is unclear if the majority of dirrs ever undergo the massive star- 
bursts that are necessary to drive such a model as photometric studies of BCDs suggest that the 
structural prope rties of these systems are fundamentally different than more “no r mal” low-surface 
brightness dirrs (jPapaderos et al.lll996l : lDoubher et al.lll999l : iMarlowe et al.lll999l : ISalzer fc Norton 
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Log(M.[Mj) 

Fig. 7.— NIR-L-Z relationship and M-Z relationship. Symbols and colors are the same as Figure 
The two fits shown are the fit to the primary comparison sample only (solid line) and fits to the 
sample including the SHIELD galaxies (dashed line). 
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I 999 I : iJanowiecki &: Salzeiil2014l ). Thus other evolutionary scenarios may need to also be consid¬ 
ered to explain the observed trends at the low-mass (and low-luminosity) end of the L-Z and M-Z 
relationships. 


As an alternative to bursty star-formation histories for isolated dirrs. iGavilan et al.l (j2013|) suggest 
that many observational features of isolated dirrs can be explained without invoking winds or large 
scale outflows, but instead rely on continuous infall or cooling of primordial gas and continuous 
star formation throughout the lifetime of the galaxy. This model provides interesting implications 
in the context of the results presented in this paper. 


In the scenario proposed bv iGavilan et al.l (j2013l i the galaxies are modeled as a single region and the 
star formation rate is controlled only by the amount of gas available and an assumed star formation 
efficiency. This model does not allow for “bursts” of rapid star formation but instead produces a 
smooth, continuous (though not constant) star-formation history for the galaxy. Their results 
suggest that with moderate to low star-formation efficiencies this continuous type of evolution 
can plausibly r eprod uce the observed dirr trends and the luminosity-metallicity relationship. The 
Gavilan et al.l (|2013l l models imply that the observed dispersion in abundance seen at a given B- 
band luminosity is a result of differing star-formation histories, owing to different star-formation 
efficiencies and/or different gas-collapse timescales. This is intriguing in the context of the work 
we have presented here, because it suggests the offset of the SHIELD sample seen in the H-band 
L-Z relationship may be indicative of different global star-formation history compared to optically 
selected samples. Since the SHIELD sample has specifically been chosen to represent some of the 
lowest-luminosity and lowest-surface brightness galaxies, it is not unreasonable to assume that our 
selection may be biased to select galaxies with different star-formation histories than those selected 
from optical catalogs. 


Tassis et al.l (120081 1 produced similar results to IGavilan et al.l (l2013li suggesting that galaxian out¬ 
flow winds are not necessary to reproduce the observed L-Z trends and attribute the inefficient 
star-formation in low-mass systems to inefficient gas cooling (i.e. longer collapse timescales) be¬ 
cause of the more extended neutral gas distributions often observed in lower mass systems. Since 
there are only a limited number of low-mass dIrr samples where high quality optical and H I data 
currently exist, it is very difficult to analyze whether this scenario is observationally consistent. 
However, we suggest this is an idea worthy of further investigation. Further studies with the 
SHIELD project, and similar projects focusing on low-mass ALFALFA sources, should lead to a 
better understanding of the nature of the gas distributions in low-mass systems and the role that 
gas dynamics play in their chemical evolution. 


The SHIELD galaxies all seem to be “offset” above the fit line from the comparison sample in 
the H-band L-Z relationship. When the “corrected” abundances are used, the offset, with respect 
to the fit from the comparison sample, is reduced, yet all but one of the SHIELD galaxies still 
appears above the fit line. This offset becomes even less pronounced when is replaced with 
Af 4.5 and disappears when stellar mass (M*) is used instead. This indicates that, on average, the 
mass-to-(i?-band) light ratio for the SHIELD galaxies differs from the comparison sample. This is 
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plausibly consistent with our previous suggestion, that our selection criteria may have biased our 
sample toward galaxies with a certain star-formation history. However, the general trend simply 
suggests that the SHIELD galaxies are too massive for their H-band luminosity. This may arise 
because at least some portion of the SHIELD galaxies do not appro ach the quasi-co ntinuous star- 
formation history that has been suggested for more massive dirrs by van Zeel (200l|). A stochastic 


star-formation history that is similar to that suggested in Ivan Zed (j200l|) but not as continuous 


(owing to the lower mass of the system) could cause the H-band luminosity to be less reflective of 
the current stellar mass in these galaxies. 


The results presented in this paper do not necessarily clarify any of the issues associated with 
the chemical evolution of dirrs, but they do suggest that varied star-formation histories need to 
be considered along with the effects of enriched gas outflows in order to adequately explain the 
observed trends. Additionally, our data suggest that while the M-Z relationship may be more 
fundamental, understanding the relationship of the mass-to-light ratio could reveal details of the 
star-formation and chemical evolution in these systems. This sample is just a first step in beginning 
to fill in the very low-luminosity range of dwarf irregulars. 


8. CONCLUSIONS 


We have used long-slit spectra from the Mayall 4m telescope at KPNO to determine the oxygen 
abundances for 8 of the 12 SHIELD galaxies. Eor two galaxies we were able to measure the 
temperature sensitive [O III] A4363 line and thus directly determine the electron temperature for 
at least on e H II regi o n. W e calculated abundances for the other 6 galaxies using the strong-line 
method of iMcGauehl (jl99l|). We calculated a n emp irical correction to the McGaugh abundance 
scale using the dirr sample from iBerg et al.l (j2012l L along with the two H II regions for which 
we detected [O III] A4363, allowing us to put our McGaugh abundances on a consistent scale 
with those that have a directly determined Tg. The galaxies in our sample have well-constrained 
distances (HST; McQuinn et ah 2014), optical luminosities (WIYN 3.5m BVR; N. Haurberg et al. in 
preparation), and NIR luminosities (Spitzer 3.6 and 4.5/im; Cannon, Marshall et al. in preparation). 
We compared the L-Z and M-Z relationships derived from our results with a substantial comparison 
set from the literature that also have well-determined abundances, dis tances, and optical and NIR 
luminosities (jLee et al.ll2006l : Ivan Zee &: Havnesll2006l : lBerg et al.ll2012l . and references therein). We 
additionally included dIr r galaxies with well-d etermined abundances that were selected from the 
ADBS blind H I-survey (jHaurberg et al.ll2013l l in our comparison sample. While more luminous 
than our data set, the latter sample is particularly interesting because the method of selection was 
very similar to that used for the SHIELD data set. 


When the L-Z relationship was examined using the R-band luminosity the SHIELD galaxies appear 
offset from the comparison sample, but that offset disappears when derived stellar mass is used 
instead. We suggest that this may indicate that a range of star-formation histories and chemical 
evolution scenarios for dwarf irregulars should be considered. The possible effects that such scenar- 
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ios may have on the the luminosity-metallicity and mass-metallicity relationships warrant further 
examination and exploration. While the M-Z relationship may be more well-constrained, the L-Z 
relationship may prove observationally important for disentangling the origins of that relationship. 
A larger and more comprehensive sample is needed and a more exhaustive comparison between 
an array of models and multiple observational parameters should be done before any model is too 
heavily favored. Such a comparison is beyond the scope of this work. 


In this work, we identified two new XMD galaxies in the local universe, AGC 112521 and AGO 
111164, and sho wed that blind H I- surveys like ALFALFA are effective in providing candidate 
XMD dirrs (e.g., ISkillman et al.ll2013l ). As follow-up studies of the low-H I mass sources from the 
ALFALFA survey continue, a larger sample of low-luminosity gas-rich dirrs will become available 
and allow us to start to unravel some of the puzzles concerning the chemical evolution of dwarf 
galaxies. Additionally, current work on resolved stellar populations in the SHIELD galaxies and 
the Leo P dwarf promise to further illuminate the star-formation histories in these very low-mass 
galaxies and lead to a deeper understanding of galaxian chemical evolution. 
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